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Abstract. The diffraction slope parameter is investigated for elastic proton-proton and proton-antiproton 
scattering based on the all available experimental data at low and intermediate momentum transfer va- 
lues. Energy dependence of the elastic diffraction slopes is approximated by various analytic functions. 
The expanded "standard" logarithmic approximations with minimum number of free parameters allow 
to describe experimental slopes in all available energy range reasonably. The estimations of asymptotic 
shrinkage parameter d v were obtained for various \t\ domains based on the all available experimental 
data. Various approximations differ from each other both in the low energy and very high energy domains. 
Predictions for diffraction slope parameter are obtained for elastic proton-proton scattering at NICA, RHIC 
and LHC energies, for proton-antiproton elastic reaction in FAIR energy domain for various approximation 
functions. 

PACS. 13.75.Cs Nucleon-nucleon interactions - 13.85.Dz Elastic scattering 
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1 Introduction 

Elastic hadron-hadron scattering, the simplest type of had- 
ronic collision process, remains one of the topical the- 
oretical problems in fundamental interaction physics at 
present. Forward elastic scattering process is an excellent 
test for some fundamental principles (unitarity, analytic- 
ity, asymptotic theorems) of modern approaches. In the 
case of pp and pp elastic scattering, although many ex- 
periments have been made over an extended range of ini- 
tial energies and momentum transfer, these reactions are 
still not well understood. One can suggests that by the 
time the accelerator complexes like RHIC, LHC etc. are 
operating, the interest to the soft physics increases sig- 
nificantly. In the absence of a pure QCD description of 
the elastic pp/pp and these large-distance scattering states 
(soft diffraction), an empirical analysis based on model- 
independent fits to the physical quantities involved plays 
a crucial role pQ. Therefore, empirical fits of energy depen- 
dencies of global scattering parameters have been used as 
an important source of the model-independent informa- 
tion. This approach for a to t and p was recently used in 
[2] . The third important quantity for nucleon elastic scat- 
tering is the slope parameter. The nuclear slope B for elas- 
tic scattering is of interest in its own right. This quantity 
defined according to the following equation: 



In 



da (s, t) 
dt 



(1) 



is determined experimentally by fitting the differential 
cross section da/dt at some collision energy yfs. On the 
other hand the study of B parameter is important, in par- 
ticular, for reconstruction procedure of full set of helicity 
amplitudes for elastic nucleon scattering [2] . In the last 20- 
30 years, high-energy pp colliders have extended the max- 
imum pp collision energy from y/s ~ 20 GeV to y/s ~ 2 
TeV, the RHIC facility allows to obtain pp data up to 
= 500 GeME As consequence, the available collection 
of pp and pp slope data from literature has extended. The 
present status of slope for elastic pp and pp scattering 
is discussed over the full energy domain. Predictions for 
further facilities are obtained based on the available ex- 
perimental data. 



2 Experimental slope energy dependence 

We have attempted to describe the energy behavior of 
the elastic nuclear slopes for pp and pp reactions. The 
classical Pomeron theory gives in first approximation the 
following expression for the differential cross section of 
elastic scattering at asymptotically high energies: 

da/dt ocs^W- 1 ), 

where a-p(t) is a Pomeron trajectory. If a-p(t) is linear 
function of momentum transfer, i.e. a-p(t) = a-p(O) + a'-pt, 
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then for the slope parameter B (s) at some t using the 
definition ([!]) one can obtain 

B (s) oc 2a'-p In s. 

Indeed the ensemble of experimental data on slopes for 
elastic nucleon collisions can be approximated reasonably 
by many phenomenological approaches, at least for y/s > 
20 GeV. But models contradict the experimental data at 
lower energies and / or phenomenological approaches have 
a significant number of free parameters as usual. On the 
other hand it is apparent from previous investigations that 
the experimental data on slope parameters do not follow 
the straight lines at any initial energies when plotted as 
function of Ins. The new "expanded" logarithmic param- 
eterizations with small number of free parameters have 
been suggested in [5][3] for description of the elastic slope 
at all available energies. Thus taking into account stan- 
dard Regge parametrization and quadratic function of log- 
arithm from [4] the following analytic equations are under 
study here: 

B(s,t) = B + 2 fll In (s/«o), (2a) 
B (s, t) = B Q + 2ai In (s/s ) + a 2 [In (s/s )P , (2b) 
B (s, t) = B + 2oi In (s/s ) + a 2 (s/s a ) a3 , (2c) 
B (s, t) = B + 2a, In (s/s ) + a 2 [In (s/s Q )} 2 , (2d) 

where So = 1 GeV 2 , in general case parameters Bo, a,i,i = 
1 — 3 depend on range of \t\ which is used for approxima- 
tion. The function (|2dj) is a special case of (2b) at fixed 
value 03 = 2. Additional terms in (]2b[) - (|2d|) take into ac- 
count the non- logarithmic part of the energy dependence 
of the elastic nuclear slopes. Approximation function (|2cp 
is analogy of parametrization of momentum slope depen- 
dence from [5]. One can see the parametrization (|2c|) is 
compatible to first order with the Regge pole model where 
the additional term represents the interference between 
the Pomeron and secondary trajectories [5]. 

Most of experimental investigations as well as theoreti- 
cal models are focused on the diffraction region \t\ ~ 0—0.5 
GeV 2 . In this paper we study all available experimental 
data for nuclear slope parameter up to \t\ ~ 3.6 GeV 2 . Ex- 
perimental data are from [6.7,8. 9 TOlfTTj . The data sam- 
ple consists of about 500 experimental points. The en- 
semble of experimental data has been slightly specified 
and improved in comparison with [3]. In particular, re- 
sults obtained in [9] for form-factor parametrization of 
da/dt were included in the data sample after detail ad- 
ditional investigation. The total number of experimental 
points is equal 142/138 for pp/pp scattering at low 
respectively. In the intermediate \t\ domain experimental 
data set is 136/69 for pp/pp reaction respectively. Thus the 
experimental sample is significantly larger than that for 
some early investigations [2,4,5,6, 12] . The careful analy- 
sis of data sample allows to suggest that the influence of 
double counting in the experimental data is negligible. It 
should be emphasized that the experimental data for in- 
termediate \t\ range were separated on two samples which 
corresponded the various parametrization types for differ- 
ential cross-section, namely, linear, hi(da/dt) oc (— B|t|), 



and quadratic, In (da/dt) oc (— B\t\ ± Ct 2 ), function. Here 
B, C > are suggested. As known the measurements of 
nuclear slope, especially at intermediate \t\ do not form 
a smooth set in energy, in contrast to the situation for 
global scattering parameters p and a to t, where there is a 
good agreement between various group data [4j. Detailed 
comparisons of slope data from different experiments are 
difficult because the various experiments cover different 
\t\ ranges, use various fitting procedures, treat systematic 
errors in different ways, and, moreover, some experimen- 
tal details are lost, especially, for very early data. We have 
tried to use as much as possible data for fitting from avail- 
able samples. But some of the B values were not further 
used, either due to internal inconsistencies in the fitting 
procedure, or as redundant in view of a better determi- 
nation at a nearby initial energy. Thus the data samples 
for approximations are somewhat smaller because of ex- 
clusion of points which, in particular, differ significantly 
from the other experimental points at close energies. Crit- 
ical for a consistent determination of slope parameters is 
the choice of the range |i| mnl < \t\ < |t|max over which 
the fit of da/dt is performed. It seems both the mean 
value of \t\ (|i|) and |i|-boundaries of corresponding mea- 
surements are important for separation of experimental 
results on different |i| domains. Here the \t\ are calculated 
with taking into account approximations of experimental 
da/dt distributions instead of identifying of |i| with mean 
point of |i|-range as previously [3]. Errors of experimental 
points include available clear indicated systematic errors 
added in quadrature to statistical ones. One need to em- 
phasize the systematic errors caused by the uncertainties 
of normalization (total or/and differential cross sections) 
are not taken into account if these uncertainties are not 
included in the systematic errors in the original papers. 

Let us describe the fitting algorithm in more detail. We 
use the fitting procedure with standard likelihood function 
for this investigation of nuclear slope parameter. In accor- 
dance with p~2] let us define the quantity 

A xi (,„ a)s ( B: »-^ M f, (3) 

where B l m is the measured value of nuclear slope at Si, 
B (s^, a) is the expected value from the one of the fitting 
functions under consideration, ai is the experimental error 
of the i-th measurement. The parameters ctj are given by 
the TV-dimensional vector a = {a±, a at}. Our fitting al- 
gorithm is some similar to the " sieve" algorithm from [12] 
with following modification. We reject the points which a 
priori differ significantly from other experimental data at 
close energies. The step allows us to get a first estimation 
of x 2 /n.d.f. with minimum number of rejected points. The 
fit quality is improved at the next steps consequently. As 
indicated above smoothness of experimental slope energy 
dependence differs significantly for data samples in vari- 
ous |i|-domains and for various parameterizations of da/dt 
(see below). The Axf (s^a) absolute value can be large 
for one data sample but it can be acceptable for another 
sample at the same time. Therefore we suggest to use the 
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relative quantity 



Ax} {sj\a) 
X 2 /n.d.f. 



(4) 



in order to reject the outliers (points far off from the fit 
curve to the data sample) instead of constant cut value 
Axf (si; oOmax fr° m the "sieve" algorithm [12]. One needs 
to emphasize the fit function with best x 2 / n -d.f. among 
(|2a|) - ([2d)l is used for expected value calculation in ([3]) 
at the each algorithm step. The points with n x > n max 
are excluded from future study in our algorithm, where 
the n™ ax is some empirical cut value. The conventional 
fit is made to the new "sifted" data sample. We consider 
the estimates of fit parameters as the final results if there 
are no excluded points for present data sample. We use 
the one value n™ ax = 2 for all data samples considered in 
this paper below. The fraction of excluded points is about 
2% for pp as well as for pp elastic scattering for low \t\ 
domain. The maximum relative amount of rejected points 
is about 3%/12% for linear \ndo/dt parametrization and 
6%/15% for quadratic one at intermediate |t| values for 
pp/pp scattering respectively. 



2.1 Low \t\ domain 

The energy dependence for experimental slopes and corre- 
sponding fits by (j2aj) - (|2d"|) are shown at FigQ] and Figj2] 
for pp and pp correspondingly. The fitting parameter val- 
ues are indicated in Table 1. One can see that the fitting 
functions (|2aj) . ([2d)l describe the pp (Figfl]) andpp (Figf2]) 
experimental data statistically acceptable only for y/s > 5 
GeV. Additional study demonstrate that the extension of 
approximation domain to the lower energies for parame- 
terizations (f2aj) and (I2djl results in significant increasing 
of x 2 /n.d.f. for pp data. Thus these fit functions allow to 
get a statistically acceptable fit qualities only at yfs > 5 
GeV for pp scattering as well as for elastic pp reaction. The 
parameter a\ for Regge-inspired function (|2a[) is close to 
estimation for Pomeron parameter a' v « 0.25 GeV -2 for 
pp data but this parameter is some larger than a' v esti- 
mate for proton-proton data. The a^ p value for (|2a| from 
Table 1 is equal within the errors the early experimental 
estimations of "true" Pomeron shrinkage parameter [13] . 
The parameter a\ for function (|2b|) is equal within errors 
to the above estimation for Pomeron theory parameter a'p 
for pp data. The value of a pp parameter for fitting function 
(12cj) agrees with early value for similar fit of slope momen- 
tum dependence [5] but the present result is more precise 
than previous one. The RHIC point for proton-proton col- 
lisions has a large error and can't discriminate the approx- 
imations. Fitting functions (|2b[) . (|2"c|) allow us to describe 
experimental data at all energies with reasonable fit qual- 
ity for pp. The functions (|2a|) - (l2c|) are close to each other 
at -^/s > 5 GeV, especially, the Regge model approxima- 
tion and (pc|) fit. It seems the ultra- high energy domain is 
suitable for separation of various parameterizations. The 
qualities of ([2b|) , (f2"c|) approximations for pp elastic scatter- 
ing data are much poorer because of very sharp behavior 
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Fig. 1. Energy dependence of the elastic slope parameters 
for proton-proton scattering for low \t\ domain. Experimental 
fitted points are indicated as •, unfitted points are indicated as 
*. The curves correspond to the fitting functions as following: 
(l2al) - dot, (f2bl) - thick solid, OB - dot-dashed, J2dl) - thin 



of experimental data near the low energy boundary. But 
one can see that the functions (|2b| , (|2c|) agree with exper- 
imental points at qualitative level and (very) close to each 
other for all energy range. The fit quality is some better 
for function (|2b[) than that for parametrization (|2"c|) both 
for pp and pp data. Additional study of antiproton-proton 
data shows that the increasing of low boundary of fit range 
(s mm ) leads to the better fit quality for functions (|2bp. 
(|2c|) but at the same time, obviously, to the lost of some 
low-energy pp data. The fit quality changes dramatically 
at small increasing of s mm from low boundary value 4m 2 
to 3.72 GeV 2 . It was obtained x 2 /n.d.f.~ 5.4/5.5 for func- 
tion (j2b|) / (|2"c|) respectively for fit range s mul > 3.72 GeV 2 . 
On the other hand the data sample is about 75% from 
maximum one in this case. Thus it seems the s m j n = 3.72 
GeV 2 one of the optimum values from point of view both 



fit quality and closing to the threshold s r 



4m 2 . The 



average value of © for excluded points is equal 15.6 for 
pp and 649.8 for pp data sample for parametrization (2b). 

One can get predictions for nuclear slope parameter 
values for some facilities based on the results shown above. 
The B values at low \t\ for different energies of FAIR, 
NIC A, RHIC, and LHC are shown in the Table 2. As ex- 
pected the functions (|2b|) and (|2c|) predicted very close 
slope parameter values for FAIR. The approximation func- 
tion (|2"a|) and (|2d[) can predict for y/s > 5 GeV only. Func- 
tions (|2ajl - (|2c]) predict much smaller values for B in high- 
energy pp collisions than (|2d[) approximation. Perhaps, the 
future more precise RHIC results will agree better with 
predictions based on experimental data fits under study. 
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Table 1. Fitting parameters for slope energy dependence in low \t\ domain 



Function 






Parameter 








B , GeV~ 2 


oi, GeV" 2 


o 2 , GeV~ 2 




X 2 /n.d.f. 


proton-proton scattering 


l2aj) 


8.41 ±0.08 


0.271 ± 0.007 






179/95 


dig 


8.73 ±0.12 


0.250 ± 0.009 


-28 ±2 


-3.50 ±0.14 


349/135 




8.31 ±0.08 


0.279 ± 0.007 


-188 ±23 


-2.24 ±0.09 


358/135 


(2d) 


9.3 ±0.3 


0.11 ±0.06 


0.03 ±0.01 




171/94 


proton-antiproton scattering 


(2a) 


10.0 ±0.2 


0.215 ±0.012 






32/27 


(f2bl> 


12.02 ± 0.06 


0.121 ±0.006 


495 ± 98 


-13.0 ±0.5 


1220/132 




12.06 ±0.05 


0.119 ±0.005 


(3.1 ±0.4) • 10 6 


-9.36 ±0.10 


1303/132 


(Hd| 


11.4 ± 1.0 


0.05 ±0.11 


0.017 ±0.011 




29/26 



Table 2. Predictions for nuclear slope based on the parameterizations pall - (|2d[) for low \t\ domain 



Fitting 
function 










Facility 


energ 


ies, y^s 












FAIR, GeV 




NICA, GeV 


RHIC, TeV 




LHC, TeV 




3 


5 


6.5 


14.7 


20 


25 


0.2 


0.5 


14 


28 


42* 


12a| 




11.38 


11.61 


12.31 


11.66 


11.90 


14.15 


15.15 


18.76 


19.51 


19.95 




12.57 


12.80 


12.93 


13.32 


11.67 


11.91 


14.02 


14.94 


18.28 


18.97 


19.37 


mi 


12.59 


12.83 


12.95 


13.34 


11.65 


11.90 


14.22 


15.25 


18.96 


19.74 


20.19 






11.90 


12.01 


12.43 


11.70 


11.96 


15.00 


16.67 


24.44 


26.39 


27.58 



The ultimate energy upgrade of LHC project [2]. 
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Our prediction with (I2djl function for ultimate energy of 
LHC agrees well with early prediction for close SSC energy 
based only on slope data in the region 5 < ^fs < 62 GeV 
[4]. It should be emphasized that the fits (|2a|) - (|2c|) of 
available experimental data predict the slower increasing 
of B with energy than most of phenomcnological models 
[15]. The B values predicted for LHC at ^/s = 14 TeV 
by (|2a| and (|2c|) are most close to the some model ex- 
pectations ,16,17]. Moreover one needs to emphasize that 
the model estimates at y/s = 14 TeV were obtained for 
B (t = 0) and the ^-dependence of slope shows the slight 
decreasing of B at growth of momentum transfer up to 
\t\ m 0.1 - 0.2 GeV 2 in particular for the models [TOUT?] . 
Thus one can expect the some better agreement between 
model estimations and predicted values of B from Table 
2 for finite (non-zero) low \t\ values. Possibly the satu- 
ration regime, Black Disk Limit, will be reached at the 
LHC. The one of the models in which such effects appear, 



namely, DDM predicts the slope B (t 
s/s = U TeV [18]. 



0) w 23 GeV -2 at 



Fig. 2. Energy dependence of the elastic slope parameters for 
antiproton-proton scattering for low \t\ domain. Experimental 
fitted points are indicated as ■. unfitted points are indicated 
as *. The correspondence of curves to the fit functions is the 
same as previously. The inner picture shows the experimental 
data and fits at the same scale as well as for fig.l 



2.2 Intermediate \t\ domain 

As indicated above the situation is more complicated for 
intermediate \t\ domain. Differential cross section is ap- 
proximated by linear, In (da /dt) cx (— 2?|£|), or / and quad- 
ratic, In (da /dt) oc (— B\t\ ± Ct 2 ), function in various ex- 
periments, \t\ ranges used for da / dt approximations differ 
significantly etc. For quadratic exponential parametriza- 
tion the B and C parameters are highly correlated by fits. 
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Fig. 3. Energy dependence of B in proton-proton (a,b) and proton-antiproton (c,d) elastic scattering for linear (a,c) and 
quadratic (b,d) exponential parametrization of differential cross section. Experimental points are indicated as close / open 
circles (squares) for pp (pp) for (a,c) / (b,d), unfitted points are indicated as The correspondence of curves to the fit 
functions is the same as previously 



Figure [3] shows the experimental data and correspond- 
ing fits for slope parameter energy dependence at interme- 
diate \t\ for pp and pp elastic scattering. The Fig|3ii and 
FigEt correspond to the linear exponential approxima- 
tion of differential cross-section for pp and pp respectively. 
Experimental data obtained at quadratic exponential fit 
of da/dt and fitting functions (|2aj) - (|2d[) are presented 
at FiglSb for pp and at Fig|3ji for pp collisions. The fit- 



ting parameter values are indicated in Table [3] for vari- 
ous interaction types and for different da/dt parameter- 
izations. Usually the fit qualities are poorer for interme- 
diate \t\ values than that for low \t\ range in pp elastic 
collisions for linear exponential parametrization of da/dt. 
The fitting functions (|2ap and (|2d|) agree with experi- 
mental points qualitatively both for linear (Fig(3ji) and 
quadratic (Fig{3]b) exponential parameterizations of da/dt 
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for y/s > 5 GeV only. The "expanded" functions J2EJ, (J2cJ) 
approximate experimental data at all energies reasonably 
with very close fit qualities. But the (|2b[) function shows a 
very slow growth of slope parameter with energy increas- 
ing at y/s > 100 GeV (Fig[3k.). It should be stressed that 
the fitting function (|2d[) predicts decreasing of the nuclear 
slope in high energy domain. Such behavior is opposite 
the other fitting function (f2"a[) - (|2"c)) . In the case of linear 
exponential approximation of da/dt mean value of char- 
acteristic (|4| is about 5.3 for excluded pp data points with 
(|2b[) function and ~ 2.6 for points excluded from pp 
fitted data sample for (|2dj) fitting function. 

One can see that the pp experimental data admit the 
approximation by (|2d[) for all energy range but not only for 
•y/i > 5 GeV. Indeed the fit quality for the first case much 
better than for second one. Additional analysis demon- 
strated just the same behavior of fit quality for function 
(|2a|) too. Thus pp experimental points from linear expo- 
nential parametrization of differential cross-section are fit- 
ted by (|2a| and (|2d[) at all energies not only for yfs > 5 
GeV. The parameter values are shown in Table [3] for ap- 
proximation by (j2"aj) , (|2d[) of all available experimental 
data. The pp data disagreement with Regge-inspired fit- 
ting function very significantly (FigJSJ;). Functions (f2"c)) 
and (j2d|) show a very close behavior at all energies for pp 
data from linear parametrization of In da/dt. These fit- 
ting functions have a better fit quality than (|2bj) but fits 
(|2c|) . (12d[) are still statistically unacceptable. As previously 
experimental data at Fig(3]d allow the approximation by 
([2a]) and (|2d|) for all energy range but not only for y/s > 5 
GeV. The fit qualities are better in the first case of energy 
range and fitting parameters are indicated in the Table 
[3] for this energy range namely. As above the functions 
(|2c|) and (|2d[) show a close fit quality which is some bet- 
ter than this parameter for (|2b[) fitting function. One can 
see the fit qualities for (f2"a|) - (|2d[) are better significantly 
for data from quadratic exponential parametrization of 
differential cross-sections than for data from linear ex- 
ponential approximation of da/dt. Thus fitting functions 
(|2"c|) . (|2d[) agree with data points at quantitative level for 
quadratic (FigJSU) parametrization of proton-antiproton 
In da/dt and these fits are statistically acceptable. Ex- 
cluded points are characterized by ~ 17.9 for pp data 
with (|2bp fitting function and by ~ 12.1 for pp data 
sample at (|2d|) function. 

From the quadratic exponential parametrization of dif- 
ferential cross-section one may compute the local slope at 
a certain |i|-value via the following relation 

b(s,t) =B±Cln\t\, B,C>0 (5) 

This characteristic can be useful for elastic scattering for 
study of da/dt in wider \t\ range. It is suggested b > 
according to the definition ([5]) . The fo-parameter is named 
slope too, it is evaluated at \t\ = 0.2 GeV 2 usually. One of 
the advantages of this characteristic is the expectation of 
more smooth energy (momentum) dependence than that 
for B-parameter discussed above. Indeed we have included 
the 100% of available experimental points in fitted sample 
for pp elastic scattering. But the number of points is some 
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Fig. 4. Energy dependence of the slope parameter b at |f| = 
0.2 GeV 2 for proton-proton scattering. The correspondence of 
curves to the fit functions is the same as previously 

smaller than that for S-parameter because of absent C- 
parameter values for some low energy measurements from 
[6j. We excluded one point at \/s = 1.8 TeV [19] from 
fitted sample for pp elastic reaction because there are un- 
acceptably large errors (relative error is 8b ~ 2.72) for this 
point. 

Experimental values of b depend on energy collisions 
and corresponding fits are shown at Fig|4] for pp elastic 
scattering and at Fig [5] for pp collisions. In the last case 
fit qualities for (f2"aj) . (|2d[) functions are better for fitting 
at \fs > 5 GeV only than that for fitting of all available 
energy domain. The fit parameter values are shown in Ta- 
ble |H Fit qualities are better significantly than that for 
corresponding fits of B parameter with the exception of 
(|2a| for pp data. Functions (|2a|. (|2d|) approximate b (^s) 
for pp data statistically acceptable for y/s > 5 GeV only. 
Functions (|2b|) and (|2c|) show acceptable close fit quali- 
ties and difference at high energies only. The shrinkage 
parameter for best approximation function (|2c[) is in 
a good agreement with a early results [5]. Function (|2bl) 
shows a best fit quality for pp data. Thus the " expanded" 
parameterizations (|2b[) and ([2c]) suppose statistically ac- 
ceptable representation of all available experimental data 
for b parameter both in pp and pp elastic reactions. 

We have obtained predictions for nuclear slope param- 
eters B and b for some facilities based on the fit results 
shown above. The predicted B values at intermediate \t\ 
are calculated on the base of fitting parameters obtained 
for linear exponential parametrization of da/dt. Slope val- 
ues are shown in the Table[5]for different energies of FAIR, 
NIC A, RHIC, and LHC. According to the fit range func- 
tion (|2aj) can predicts the B value for pp scattering at all 
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Table 3. Fitting parameters for slope energy dependence at intermediate \t\ 




Function 


Parameter 






B , GeV~ 2 ai, GeV~ 2 a 2 , GeV~ 2 a 3 


X 2 /n.d.f. 


proton-proton scattering, experimental data for da/dt = Aexp (—B\t\) 


l2a|) 

dig 

(2d) 


7.95 ±0A2 0.184 ±0.009 

9.7 ±0.3 0.08 ± 0.02 -22.1 ±1.4 -2.34 ±0.12 
8.51 ±0.14 0.144 ±0.010 -71 ±5 -1.49 ±0.06 

4.9 ±0.6 0.73 ±0.10 -0.09 ± 0.02 


110/34 
240/60 
240/60 
81/33 


proton-proton scattering, experimental data for da/dt — Aexp (— B \t\ ± Ct A ) 


((2aJ) 

(EH 


7.4 ±0.2 0.31 ±0.03 

9.6 ±2.4 0.16 ±0.13 -7.2 ±5.4 -1.5 ±1.0 
7.9 ±0.5 0.26 ±0.05 -23 ±16 -1.5 ±0.5 
4.1 ±0.9 1.0 ±0.2 -0.15 ±0.04 


115/33 
227/62 
228/62 
102/32 


proton-antiproton scattering, experimental data for da/dt = Aexp (-B t ) 


l2a|) 

(Eg 

(E3> 


11.25 ± 0.06 0.126 ±0.004 
(-5.0 ±0.9) ■ 10 3 0.64 ±0.02 (5.1 ± 0.9) ■ 10 3 (-1.5 ± 0.3) • 10~ 3 
-893 ±69 5.9 ±0.3 908 ± 69 (-1.44 ± 0.06) • 10~ 2 
15.46 ±0.15 -0.59 ±0.02 0.083 ± 0.003 


1111/41 

355/39 
199/39 
197/40 


proton-antiproton scattering, experimental data for da/dt = Aexp (— B\t\ ± Ci 2 ) 


Urn 
(Eg 
(H3 


11.1 ±0.2 0.171 ±0.011 
(-1±3)-10 3 0.26 ±0.06 (1±3)-10 3 (-1±4)-10~ 3 
(-2.2 ± 1.8) • 10 2 1.7 ±0.8 (2.3 ± 1.8) • 10 2 -0.015 ± 0.006 
12.7 ±0.8 -0.05 ±0.11 0.023 ±0.011 


17.3/15 
14.5/13 
13.2/13 
13.2/14 



Table 4. Fitting parameters for b energy dependence 



Function 






Parameter 








Bo, GeV~ 2 


ax, GeV~ 2 


a 2 , GeV~ 2 




X 2 /n.d.f. 


proton-proton scattering 


(2a) 


6.7 ±0.2 


0.32 ± 0.02 






25.6/25 


dig 


8.4 ± 0.4 


0.20 ±0.03 


-10.4 ± 1.4 


-2.0 ±4.5 


41.5/33 




7.6 ±0.8 


0.25 ±0.07 


-20 ± 15 


-1.2 ±0.4 


41.4/33 


(f2d|) 


5.8 ± 0.9 


0.5 ±0.2 


-0.05 ±0.04 




24.3/24 


proton-antiproton scattering 


<l2a|) 


9.2 ±0.2 


0.22 ±0.02 






40.5/13 


dig 


-234 ± 143 


0.71 ±0.07 


252 ± 152 


-0.04 ±0.02 


8.8/15 




2.4 ± 1.6 


0.46 ±0.11 


18.5 ±2.3 


-0.29 ±0.11 


9.0/15 


(Eg 


14 ± 1 


-0.47 ±0.12 


0.07 ±0.01 




8.1/12 



Table 5. Predictions for slope parameters based on the functions (|2ap - (|2d[) for intermediate \t\ domain 



Fitting- 










Facility energies, yfs 










function 




FAIR, GeV 




NICA, GeV 


RHIC, TeV 




LHC, TeV 






3 


5 


6.5 


14.7 


20 


25 


0.2 


0.5 


14 


28 


42* 












B-paramet 


er 










(2a) 


11.80 


12.06 


12.19 


12.60 


10.15 


10.32 


11.85 


12.52 


14.98 


15.45 


15.79 


(2b) 


14.08 


12.56 


12.12 


11.52 


10.31 


10.43 


11.26 


11.57 


12.64 


12.86 


12.98 


(2c) 


13.29 


12.52 


12.19 


11.50 


10.23 


10.36 


11.56 


12.09 


14.01 


14.41 


14.64 


(2d) 


13.26 


12.51 


12.19 


11.50 


10.45 


10.61 


10.38 


9.30 


0.33 


-2.53 


-4.36 


6-parameter 


(2a) 


10.00 


10.46 


10.70 


11.43 


10.52 


10.81 


13.44 


14.60 


18.81 


19.68 


20.20 


(2b) 


14.49 


12.60 


12.03 


11.24 


10.50 


10.72 


12.55 


13.32 


16.05 


16.61 


16.94 


(2c) 


14.22 


12.65 


12.10 


11.22 


10.49 


10.72 


12.79 


13.69 


16.99 


17.67 


18.08 


(2d) 


12.56 


11.99 


11.76 


11.30 


10.50 


10.72 


11.87 


11.87 


9.27 


8.22 


7.52 



*The ultimate energy upgrade of LHC project [T4] . 
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1 10 10 2 10 3 

VS, GeV 

Fig. 5. Energy dependence of the slope parameter b at \t\ = 0.2 
GeV 2 for antiproton-proton scattering. The correspondence of 
curves to the fit functions is the same as previously 



energies under study not in yfs > 5 GeV domain only. 
As expected the functions (|2"c)) and (|2d[) predicted very 
close slope parameter values for FAIR. All fitting func- 
tions, especially |2b|) and (|2c|) . predict the close values for 
nuclear slope in NICA energy domain. Functions (|2"a|) 
(|2c|) predict larger values for B in high-energy pp colli- 
sions than (|2dj) approximation. Perhaps, the future more 
precise RHIC results will be useful for discrimination of 
fitting functions under study for intermediate \t\ values. 
The function (|2dj) with obtained parameters predicts neg- 
ative B values at future LHC energies. It should be empha- 
sized that various phenomenological models predict a very 
sharp decreasing of nuclear slope in the range \t\ ~ 0.3— 0.5 
GeV 2 at LHC energy s/s = 14 TeV [11]. Just the negative 
B value predicted for LHC at sfs = 14 TeV by (JM]) is most 
close to the some model expectations 17,20 . Taking into 
account predictions in Table [2] based on the fitting func- 
tions (|2aj) - (|2dj) for low \t\ one can suggest that the model 
with hadronic amplitude corresponding to the exchange of 
three pomerons [15p7] describes the nuclear slope some 
closer to the experimentally inspired values at LHC energy 
both at low and intermediate \t\ than other models. One 
can see the functions (|2bj) and (f2c|) predict very close val- 
ues of b up to LHC energies. Function ([2d)) shows a much 
slower decreasing of b at LHC energy domain than that for 
B parameter. The values of b parameter differ significantly 
from each other for various fitting functions in ultra-high 
energy domain for pp collisions and at low energies for 
pp elastic scattering. It seems b parameter is more per- 
spective for distinguishing of Regge-inspired function (|2"af 
from "expanded" parameterizations (|2bj) . (|2c|) at y/s > 0.5 



TeV than B because of larger differences between predic- 
tions of b for corresponding functions. 



2.3 AB and NN data analysis 

Phenomenological models predicts the zero difference of 
slopes (AB) for proton-antiproton and proton-proton elas- 
tic scattering at asymptotic energies. Here the difference 
AB is calculated for each function (|2a[) - {2d} under study 
with parameters corresponded pp and pp fits: ABi (s) = 
Bf p (s) - Bf (s) , i = 2a,...2rfl It should be stressed that 
the equal energy domain are used in pp and pp fits for 
AB calculations, i.e. the parameter values obtained by 
(|2dj) fitting function for pp data from linear exponential 
fit of da/dt for y/s > 5 GeV are used for corresponding 
AB definition. The energy dependence of AB is shown 
at FigJSk and Fig[Sh) for low and intermediate \t\ respec- 
tively. One can see that the difference of slopes decreas- 
ing with increasing of energy for low |t| domain (FigGjH). 
At present the proton-proton experimental data at high- 
est available energy 200 GeV don't contradict with fast 
(square of logarithm of energy) increasing of slope at high 
energies in general case. Such behavior could be agreed 
with the asymptotic growth of total cross section. But on 
the other hand the quadratic logarithmic function (|2dj) 
leads to very significant difference AB for pp and pp scat- 
tering in high energy domain for both low (Fig[6H) and in- 
termediate (FigEb) values of |t|. The only Regge-inspired 
function (f2a|) predicts the decreasing of AB with energy 
growth at intermediate \t\ (FiglSb). The parameterizations 
(|2b|) - ([2d|) predict the decreasing of difference of slopes at 
low and intermediate energies and fast increasing of AB 
at higher energies for intermediate \t\ domain (FigHb). As 
expected the most slow changing of AB is predicted by 
Regge-inspired function (|2a|) at asymptotic energies. All 
fitting functions with experimentally inspired parameters 
don't predict the constant zero values of AB at high en- 
ergies. But it should be emphasized that only separate 
fits were made for experimental data for pp and pp elas- 
tic reactions above. These results indicate on the impor- 
tance of investigations at ultra-high energies both pp and 
pp elastic scattering for many fundamental questions and 
predictions related to the general asymptotic properties of 
hadronic physics. 

Also we have analyzed general data samples for pp 
and pp elastic scattering. Slope parameters (B and b) 
shows a different energy dependence at y/s < 5 GeV in 
proton-proton and antiproton-proton clastic reactions in 
any \t\ domains under study. Thus slopes for nucleon- 
nucleon data are investigated only for y/s > 5 GeV below. 
We have included in fitted samples only pp and pp points 
which were included in corresponding final data samples 
at separate study pp and pp elastic reactions above. We 
did not exclude any points from NN data sample, we 
change only the low energy boundary for fitted domain. 

2 Obviously, one can suggest various combinations of fitting 
functions for AB calculations, for example, the difference be- 
tween fitting functions with best fit qualities etc. 
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Fig. 6. The energy dependence of the difference of elastic slopes for proton-antiproton and proton-proton scattering in low |t| 
domain (a) and in intermediate \t\ range for linear exponential fit of differential cross-section (b). The correspondence of curves 
to the fit functions is the same as above 



FigUJ shows the experimental data for slope in nucleon- 
nucleon elastic scattering against collision energy at low 
\t\. As seen from Fig{7] there is no experimental data for 
pp between y/s — 5 GeV and y/s — 10 GeV. This en- 
ergy domain will available for further FAIR facility. We 
have fitted the general nucleon-nucleon data sample at 
range of low energy boundary s mm = 25 - 400 GeV . 
The fitting parameter values are indicated in Table [5] on 
the first line for s mm = 25 GeV 2 and on the second 
one for s mm = 400 GeV 2 . The fit quality improves for 
all parameterizations under consideration at increasing 
of s mm . One need to emphasize the fit quality is some 
poorer (x 2 /n.d.f.~2.3 - 2.4) at yfs > 10 GeV than that 
for y/s > 5 GeV for functions (f2"aj) and (|2d[) . In the first 
case the ai parameter value is close to the Regge model 
prediction within errors for function (|2d[) only. The value 
ai = 0.24±0.01 GeV" 2 obtained for fitting functions ([2b)l 
and (f2c|) at s mm = 400 GeV 2 agrees very well with esti- 
mate for asymptotic shrinkage parameter a'-p . All func- 
tions (|2a[) - (|2d[) are close to each other at energies up 
to y/s ~ 10 TeV at least and shows quasi-linear behav- 
ior for parameter values obtained by fits with s mm = 25 
GeV 2 . Fitting functions (2a|) - j2d} are shown at FigUJ 
for s mm = 400 GeV 2 . The function ([2df decreases at 
ultra-high energies y/s > 20 TeV and clear separation of 
various fitting functions is accessible in the LHC energy 
domain. We have analyzed the nucleon-nucleon data for 
slope parameters B and b at intermediate \t\ values for 
linear and quadratic exponential parametrization of da/dt 
respectively. Fit results are shown in Table [71 Experimen- 
tal pp and pp data for B differ significantly up to y/s ~ 10 
GeV at least that results in unacceptable fit qualities for 



CM 
■ 




gU 

1 10 1 2 1 3 1 4 1 5 

Vs, GeV 

Fig. 7. Energy dependence of the elastic slope parameter for 
nucleon-nucleon scattering for low |i| domain. Experimental 
fitted points are indicated as • for pp and as ■ for pp. Fits are 
shown for ^s mm = 20 GeV. The correspondence of curves to 
the fit functions is the same as previously 
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all functions under study (x 2 /n.d.f. ~ 30.1 for best fit 
by quadratic logarithmic function). Additional analysis 
demonstrate the improving of fit quality for (|2b[) - (|2d[) 
with increasing of low energy boundary from s mm = 25 
GeV 2 up to s min = 400 GeV 2 . The values of fit param- 
eters for the last case are shown in Table [7] The Regge- 
inspired function (|2a|l contradicts with experimental data. 
We would like to emphasize that the best fit quality for 
l[2a|l is obtained at s min = 100 GeV 2 (x 2 /n.d.f. ~7.8) but 
it is statistically unacceptable too. Functions j2b| - (|2dj) 
represent experimental B (y/s) reasonably and have a very 
close fit qualities. One can see from Table [7] the statisti- 
cally acceptable fits have been obtained for b parameter at 
s min = 400 GeV 2 only. Experimental data and fit func- 
tions are presented at the Fig[8] Functions (|2bj) - (|2d|) 
show close fit qualities. Best fit is (|2dp but "expanded" 
parameterizations agree with data too. One needs to em- 
phasize the significant errors and absence of experimental 
points at y/s ~ 0.1 — 2 TeV that trouble the more clear 
conclusion. The RHIC as well as LHC data for nucleon- 
nucleon differential cross-section at intermediate \t\ will 
be helpful for distinguishing of various fit functions. 

One can conclude the slope parameters for pp and pp 
elastic scattering show universal behavior at y/s > 20 GeV 
and " expanded" functions represent the energy dependen- 
cies for both low and intermediate \t\ ranges for this en- 
ergy domain. Thus quantitative analysis of slopes at dif- 
ferent \t\ allows us to get the following estimation of low 
energy boundary: ^/s ~ 20 GeV for universality of elas- 
tic nucleon-nucleon scattering. This estimates agrees with 
results for differential cross-sections of pp and pp elastic 
reactions based on the crossing-symmetry and derivative 
relations [2]. 



3 Conclusions 

The main results of this paper are the following. Slope 
energy dependencies are analyzed quantitatively for elas- 
tic nucleon-nucleon scattering in various \t\ domains. Most 
of all available experimental data for slope parameter in 
elastic nucleon collisions are approximated by different an- 
alytic functions. 

The suggested new parameterizations allow us to de- 
scribe experimental nuclear slope at all available energies 
in low \t\ domain for pp quite reasonably. The new ap- 
proximations agree with experimental pp data at qualita- 
tive level but these fits are still statistically unacceptable 
because of very sharp behavior of B near the low energy 
limit. The best fit quality is obtained for "expanded" log- 
arithmic function both for pp and pp data. The obtained 
values of asymptotic shrinkage parameter a' v for pp elas- 
tic scattering are larger than a' v values for elastic pp re- 
actions for the same fitting functions. Various approxima- 
tions differ from each other both in the low energy and 
very high energy domains. Predictions for slope param- 
eter are obtained for clastic proton-proton and proton- 
antiproton scattering in energy domains of some facilities 
at low momentum transfer. Our predictions based on the 
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Fig. 8. Energy dependence of the elastic slope parameter 6 at 
\t\ = 0.2 GeV 2 for nucleon-nucleon scattering for intermediate 
\t\ domain. Experimental fitted points are indicated as A for 
pp and as T for pp. Fits are shown for ^s mm = 20 GeV. The 
correspondence of curves to the fit functions is the same as 
previously 



all available experimental data don't contradict the phe- 
nomenological model estimations qualitatively. The situa- 
tion is more unclear at intermediate \t\ values than for low 
\t\ domain. Only the qualitative agreement is observed be- 
tween approximations and experimental points both for pp 
and pp collisions for linear exponential parametrization of 
da/dt because of poorer quality of data. The "expanded" 
logarithmic function describes of pp data for B parameter 
for any differential cross-section parametrization reason- 
ably. Best fit quality is obtained for quadratic function 
of logarithm for pp data. One needs to emphasize that 
this function allows us to describe pp data at all avail- 
able energies and shows a statistically acceptable fit qual- 
ity for data sample obtained from quadratic exponential 
parametrization of da/dt. Slope parameter b calculated at 
\t\ = 0.2 GeV 2 shows more smooth energy dependence. We 
have obtained an acceptable fit qualities for "expanded" 
functions both for pp and pp data at all initial energies. 
The obtained values of asymptotic shrinkage parameter 
a' v for pp elastic scattering are close to the early results 
both in low and intermediate \t\ domain. As well as for low 
\t\ domain predictions for slope parameters B and b are 
obtained for elastic proton-proton and proton-antiproton 
scattering in energy domains of some facilities. It seems 
the phcnomcnological model with hadronic amplitude cor- 
responding to the exchange of three pomerons describes 
the nuclear slope some closer to the experimental fit in- 
spired values at LHC energy both at low and intermediate 
\t\ than other models. 
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Table 6. Fitting parameters for slope energy dependence in low \t\ domain for NN elastic scattering 



Function 






Parameter 








Bo, GeW 2 


ax, GeV" 2 


a 2 , GeV -2 




X 2 /n.d.f. 


(|2a[> 


8.35 ±0.06 


0.278 ± 0.005 






272/124 




7.89 ±0.11 


0.306 ± 0.008 






82.8/60 


dig 


8.21 ±0.07 


0.285 ± 0.005 


1.6 ±0.4 


-2.0 ±2.4 


271/122 




9.54 ±0.25 


0.24 ±0.01 


-514 ±244 


-3.5 ±0.3 


59.1/58 


dig 


8.05 ±0.55 


0.29 ±0.02 


0.64 ±0.55 


-0.3 ±0.4 


270/122 




9.3 ±0.2 


0.24 ±0.01 


-43 ± 18 


-0.67 ±0.08 


59.1/58 




8.48 ±0.14 


0.26 ±0.02 


(3 ±3) ■ 10 -3 




271/123 




5.6 ±0.5 


0.58 ±0.06 


-0.029 ± 0.006 




61.1/59 



Table 7. Fitting parameters for energy dependence of slope parameters at intermediate \t\ for NN elastic scattering 



Function Parameter 





B , GeV" 2 


ax, GeV" 2 


a 2 , GeV" 2 


03 


X 2 /n.d.f. 


B-parameter 


l2a| 


5.44 ±0.12 


0.352 ±0.007 






268/31 




-242 ± 43 


1.00 ±0.05 


267 ±45 


-0.06 ±0.01 


110/31 


Hi 


-135.9 ±0.5 


2.31 ±0.06 


151.0 ±0.7 


-0.040 ± 0.002 


102/31 




14.1 ± 0.7 


-0.52 ±0.07 


0.081 ±0.006 




101/32 


6-parameter 


(|2a|> 


6.97 ±0.13 


0.318 ±0.012 






221/40 




6.8 ±0.5 


0.31 ±0.03 






15.3/7 


((2b| 


7.5 ±0.2 


0.282 ±0.017 


-595 ± 662 


-6 ±3 


206/38 




-23 ±8 


0.98 ±0.19 


144 ± 40 


-1 ± 2 


2.43/5 




7.5 ±0.2 


0.281 ±0.016 


-152 ±287 


-1.7±0.6 


205/38 




9.0 ±28.5 


-1.7 ± 1.4 


13 ±31 


0.10 ±0.08 


2.40/5 


(E3) 


6.5 ±0.3 


0.40 ± 0.04 


-0.010 ±0.005 




217/39 




26 ±5 


-1.6 ±0.5 


0.17 ±0.05 




2.40/6 



The energy dependence of difference of slopes (AB) 
for proton-antiproton and proton-proton elastic scatter- 
ing was obtained for fitting functions under study. The 
AB parameter shows the opposite behaviors at high en- 
ergies for low and intermediate \t\ domains (decreasing / 
increasing, respectively) for all fitting functions with the 
exception of Regge-inspired one. The last function predicts 
the slow decreasing of AB with energy growth. It should 
be emphasized that all underlying empirical fitting func- 
tions with experimentally inspired parameter values don't 
predict the zero difference of slopes for proton-antiproton 
and proton-proton elastic scattering both at low and in- 
termediate \t\ for high energy domain. We have analyzed 
general nucleon-nucleon data samples for slopes at yfs > 5 
GeV. The "expanded" functions show the best and statis- 
tically acceptable fit qualities at yfs > 20 GeV for low \t\ 
domain. Slop analysis allows us to find the following value 
0.24 ± 0.01 for aL parameter. The estimation of asymp- 
totic shrinkage parameter a' v obtained with "expanded" 
functions for NN data agree very well with Pomeron the- 
ory expectation. The quadratic logarithmic function rep- 
resents experimental NN data for B and b slope param- 
eters with best quality. But the functions (|2b[) and (|2"c)) 
show a close qualities and agree with data reasonably. 
Therefore suggested " expanded" functions can be used as 
a reliable fits for wide range of momentum transfer at all 



energies. The universal behavior was found for available 
experimental pp and pp slopes at ^fs > 20 GeV both for 
low and intermediate \t\ that is in agreement with the hy- 
pothesis of a universal shrinkage of the hadronic diffrac- 
tion cone at high energies. 
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